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Ultra-high-frequency piecewise-linear chaos using delayed feedback loops 
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We report on an ultra-high- frequency (> 1 GHz), piecewise-linear chaotic system designed from low-cost, com- 
mercially available electronic components. The system is composed of two electronic time-delayed feedback 
loops: A primary analog loop with a variable gain that produces multi-mode oscillations centered around 
2 GHz and a secondary loop that switches the variable gain between two different values by means of a 
digital-like signal. We demonstrate experimentally and numerically that such an approach allows for the 
simultaneous generation of analog and digital chaos, where the digital chaos can be used to partition the sys- 
tem's attractor, forming the foundation for a symbolic dynamics with potential applications in noise-resilient 
communications and radar. 
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Recently, Corron et al. realized a piecewise-linear 
electronic circuit that simultaneously generates 
an analog chaotic waveform and an associated dig- 
ital waveform that represents the system's sym- 
bolic dynamics. Their technique is data-efficient 
for chaos-based applications in that the symbolic 
dynamics is sufficient for specifying the underly- 
ing deterministic features of the chaotic carrier 
signal, but their implementation was limited to 
the audio frequency range (20 Hz - 20 kHz). Here, 
we propose a new circuit architecture that oper- 
ates at ultra-high- frequencies (0.3 - 3 GHz). Our 
approach is based on two time-delayed electronic 
feedback loops, one analog and the other digital. 
Together, they ensure an increase in complexity 
of the analog signal and enable fast control over 
its dynamics, so that the overall system still be- 
haves in a piecewise-linear fashion. Our architec- 
ture represents a potential breakthrough in the 
design of inexpensive and robust electric circuits 
that generate piecewise-linear chaos. 



I. INTRODUCTION 

A picccwisc-linear (PWL) operator with two-states is 
defined as 



/(*) 



fi(x), ifx<I 
Mx), i£x>I 



(1) 



where fi(x) and f2(x) are linear operators with respect 
to x and I is a threshold. Everytime the value of x crosses 
/, f(x) switches between fi(x) and faix)- This definition 
extends naturally to multiple states and thresholds. 

Piecewise-linear operators have been used with discrete 
maps,— continuous systems,—^ and to describe certain 
classes of electronic circuits^ and regulatory systems^ 
They introduce a nonlinear component in the dynamics 
of systems (called PWL systems), a minimum require- 
ment for potentially leading to chaos. For example, when 



a PWL operator controls the parameters of discrete maps 
and linear differential equations, chaos can exist in the 
temporal evolution of the systems' state variables^ Re- 
cent studies also demonstrate that a two-state PWL op- 
erator can generate low-dimensional chaos if the system 
has hysteresis^ or uses multiple thresholds^ Such theo- 
retical evolutions, refered to as PWL chaos, can also be 
observed in experimental systems. 

Creating an experimental cmbodicment of a PWL sys- 
tem resides in the ability to design multiple linear systems 
(or a single system with variable parameters) and a phys- 
ical switch. Examples of linear systems used to produce 
experimental PWL dynamics arc LRC (inductance, resis- 
tance, capacitance) audio- frequency circuits^ They are 
popular because their parameters can be easily switched 
with conditions based on their voltage or current. This 
allows the LRC circuits to operate in two different linear 
statesAi Physical switches implemented in PWL systems 
use various components such as digital logic gates and 
operational amplifiers^ Though the switching times be- 
tween the linear states are not infinitely fast, they are 
usually considered negligible for LRC circuits (being or- 
ders of magnitude faster than the system's dynamics) and 
thus make LRC systems good for PWL experiments^^ 

A novel LRC-based implementation of a PWL sys- 
tem has been designed recently by Corron et al.r^ In 
their circuit, negative damping (gain greater than one) 
induces oscillations about a piecewise-constant voltage 
with a growing amplitude that is bounded by a piecewise 
constant driving voltage. The driving voltage switches 
its sign in response to a guard condition on the system's 
present dynamics, and leads to the generation of audio- 
frequency PWL chaos. Such chaos is fully characterized 
by the temporal evolution of the piecewise-constant volt- 
age that serves as an easily-accessible, real-time symbolic 
dynamics. This approach, with its simultaneous gener- 
ation of analog and digital chaos, demonstrates advan- 
tages over traditional chaos generators from an applica- 
tion point-of-view. One of them is the existence of a 
matched filter for robust chaos-based communications^ 
and radars^ in noisy environments. 
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DIGITAL CONTROL LOOP 




ANALOG FEEDBACK LOOP 



FIG. 1. Piecewise-linear chaotic waveform generator. An 
analog feedback loop connects the output voltage («out) of 
a variable-gain amplifier (VGA, Hittite HMC287MS8) and 
its input voltage (vi n ). The VGA output also feeds a digital 
control loop with a logarithmic amplifier (LGA, Analog De- 
vices AD8319) and a transistor-transistor-logic (TTL, Texas 
Instruments SN74AUCIG04) inverter gate in series (a second 
TTL gate is driven in parallel in order to monitor s(t) with 
correct impedance matching). We monitor the signals v(t) 
and s(t) using an oscilloscope (osc). 



The range of applications for PWL chaos, as envisioned 
by Corron et at, motivates an ultra-high-speed version 
of their circuit. However, non-ideal behaviors of LRC 
circuits make it difficult to realize an ultra-high-speed 
PWL mode-of-opcration. Above a certain frequency, the 
propagation times of signals through LRC circuits are no 
longer negligible^ and the dynamics evolve on a time 
scale that is comparable to the switching times of the 
system's electronic PWL operators. Furthermore, non- 
linearities tends to arise in LRC systems at ultra-high- 
frequency (UHF) 17 ' 18 and thus introduce undesired ef- 
fects into the PWL system. Therefore, in order to over- 
come these obstacles, a new architecture must addresses 
these challenges that are at the forefront of nonlinear dy- 
namics and electronic design. 

In this paper, we present the design and operating 
characteristics of a novel chaotic system, shown in Fig. [TJ 
that uses feedback loops to realize UHF low-dimensional 
chaos. In the following sections, we first describe the two 
main components of our electronic system: (i) an ana- 
log time-delayed feedback-loop, which generates a multi- 
mode oscillations; and (ii) a digital time-delayed con- 
trol loop, which bounds the amplitude of the oscillations 
using a threshold condition and generates a digital-like 
waveform that follows the threshold crossings. Then, we 
show that the dynamics produced are chaotic and PWL 
using a time-delayed return map. Finally, we present 
a PWL physical model that captures the dynamics dis- 
played by our chaos generator. 



II. ANALOG MULTI-MODE FEEDBACK LOOP 

The time-delayed feedback loop in See Fig. [T] is used 
to generate growing-amplitude oscillations, similar to the 
negative damping in Corron et al.'s designs . 12 ' 19 It com- 
prises a variable gain amplifier (VGA) that is self fedback. 
The VGA produces a voltage v out = v(t) that is delayed 
by time Tf by means of coaxial cables and connected to 
its input [(v ln = v(t — Tf))]. The net gain g{t) of the 
feedback loop depends on the the VGA's gain, which is 
controlled by the voltage v c t\ applied to its control port. 
If v c ti > 0.5 V then g(t) < 1, otherwise g(t) > 1. As we 
will show, these two gain values yield different dynamical 
regimes for the entire feedback loop. 

The presence of a time-delayed feedback loop has im- 
portant implications for the spectral properties of our 
system compared to that of Corron et ai. For example, it 
generates multi-mode oscillations instead of mono-mode 
oscillations.^^ Within the bandwidth B of the VGA, 
resonances that lead to multi-mode oscillations occur ap- 
proximately at integer multiples of 1 /Tf . If the pass-band 
of B is larger than the (approximate) inter-mode spacing, 
the system becomes multi-mode. 
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FIG. 2. Multimode characteristics of the feedback loop, (a) 
Experimental setup to determine the open-loop transfer func- 
tion of the feedback loop. The input of a VGA is fed by a 
sinusoidal waveform generator and a delayed version of its 
output v(t — Tf), where the signal v(t) is monitored using a 
power spectrum analyzer (PSA). For this measurement, an 
attenuator of gain g a in the feedback loop to maintain the 
net-gain g&g(t) below unity, avoiding self-oscillation and sat- 
uration at the resonant frequencies. In our experiment, we 
observe multimode characteristics with (b) v ct \ = Vl, the 
high-gain state, and (c) v ct i = «h, the low-gain state, where 
both states allow for the existence of up to ~20 modes (mode 
spacing ~l/rf ~ 25 MHz). 
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FIG. 3. Growth, saturation, and decay of v(t). (a) The feed- 
back loop is initialized with #h = 1.5 at t ~ 4/is. Electrical 
noise is amplified and band-pass filtered such that multimode 
oscillations emerge. Saturation occurs at t ~ 4.85 fj,s and then 
the gain is switched low (g^ = 0.2) at t ~ 6.25 ^s, where os- 
cillations decay in steps of rf (given by a horizontal bar in 
(d)). Zooms of v(t) showing regions where its amplitude is 

(b) governed by the system's multi-mode transfer function, 

(c) saturates, and (d) decays. Note that we choose gu and g^ 
such that the rate of growth is slower than the rate of decay 
(sh - K 1 - gi). 



In addition, the multi-mode resonances and B arc de- 
pendent on the value of v ct \. To characterize these spec- 
tral changes, we measure the feedback-loop transfer func- 
tion for two different values of the voltage control port 
w ct i, where the setup is shown in Fig. [2^,. It comprises a 
sinusoidal waveform generator with automatic frequency 
sweep, which drives the feedback loop, and a power spec- 
trum analyzer (PSA) that records the output power. 
The experiment is first realized with u ct i = «l < 0.5 V, 
shown in Fig. &p. It confirms the existence of multiple 
modes within B ~ 0.4 GHz centered at f c ~ 2.1GHz. 
For v c ti = vh > 0.5 V (Fig. [3b), we observe reso- 
nances with relatively constant amplitude, and we ap- 
proximate B ~ 0.6 GHz with a shift in the central fre- 
quency f c ~ 2.5 GHz (See Appendix). 

When the net-gain of the feedback is greater than unity 
(without an external source), the VGA amplifies multi- 
mode signals after each roundtrip in the feedback loop 
until it saturates. Figure [3] illustrates how the feedback 
loop amplifies and band-pass filters electrical noise af- 
ter propagating several time through the VGA with loop 
net-gain g{t) = gn > 1. The output voltage v(t) starts 
showing modulated oscillations (see Fig. |3p). When os- 
cillations grows beyond ^0.3 V, the VGA saturates and 
behaves nonlincarly such that the dynamics become a 
stable periodic oscillation of constant amplitude at fre- 
quency ~/ c (see Fig. |3b). We maintain the saturated 



behavior for approximately 1 //s (~25 x rf) to confirm it 
is not transient dynamics. Then, the net-gain net gain 
is switched to g(t) = gh < 1, and the oscillations in the 
feedback loop decay until they reach the noise floor of 
the VGA. 

This experiment shows that a single feedback-loop sys- 
tem can generate multi-mode UHF oscillations. Away 
from saturation, the feedback loop operates linearly and 
is governed solely by the values of its net gain (<?l,h), 
the time delay tj, and the band-pass characteristics ap- 
propriate for a given value of v ct i- This experiment also 
shows that, without proper control over the time spent in 
the growth regime, the system saturates rapidly (< 1/xs), 
which introduces nonlincarity into the dynamics. 

However, if the VGA gain is switched fast enough and 
with proper timing, it can avoid saturation and remain in 
its linear mode-of-operation. For this, in the upcoming 
section, we introduce an additional time-delayed feed- 
back loop that bounds the amplitude of the multi-mode 
oscillations. This digital control loop is specifically de- 
signed to parallel aspects of the guard condition from 
Ref. [TH, where it uses a switching conditio n 20 ! 21 that 
produces a digital- like waveform s(t). 



III. DIGITAL CONTROL LOOP 

We design our control scheme to vary v c t\ in response 
to a signal that monitors the amplitude of the multi-mode 
oscillations in v(t). When the amplitude or envelope of 
v(t) grows beyond (decays below) a set threshold, the 
control loop changes the value of v ct i and switches the 
system to its decay (growth) regime. In our design, we 
use a logarithmic amplifier (LGA) to detect the envelope 
of v(t) and a transitstor-transistor logic (TTL) inverter 
gate as a digital switch (Sec Appendix for characteriza- 
tion of these components). As a whole, the control- loop 
generates a digital signal and is designed to prevent sat- 
uration in the feedback loop. 

As shown in Fig. [TJ the feedback and control loops 
are combined by means of a voltage divider, so that v(t) 
drives the VGA, the LGA, and the TTL gate. The output 
of the TTL gate generates a digital signal s(t) that is fed 
into a fixed voltage attenuator connected directly to v c a- 
We use a fixed voltage attenuator to map the low and 
high state from s(t) (s(t) = V, and 2 V, respectively) 
onto high and low values desired for g(t) (w c ti = and 
«h, respectively). This design is adequate to realize a 
switching condition for PWL dynamics, except that it 
requires additional time-delay to guarantee proper con- 
trol on the fast dynamics existing in the feedback loop. 

With low-speed systems, the processing time of the 
components in a control loop, called the control-latency- 
time (CLT), is usually orders of magnitude faster than 
the timescales of the dynamics to be controlled and there- 
fore can be neglected^ ' 12 ! 22 In our UHF system, however, 
this is not the case. The average period in v(t) is 0.5 ns, 
which should be compared to the summed CLT of the 
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LGA and TTL gate of approximately 9 ns. As a result, 
it is not possible for the switching condition to affect in- 
stantaneously the feedback-loop oscillations as they are 
generated. 

However, we can control the envelope of the oscilla- 
tions (evolving on a slower time-scale) if s(t) is applied 
to the VGA at the proper time with respect to feedback- 
loop delay time. This is obtained by adjusting the time 
delays of both the control and feedback loops to be ap- 
proximately equal. Such an approach has been already 
used successfully to stabilize a very-high-frequency opto- 
electronic chaotic system^ 

The time delay Tf of the feedback loop is tunable with 
respect to the length of coaxial cable used in the design of 
our circuit. First, we guarantee that Tf is larger than the 
CLT and choose Tf ~ 41 ns. Then, we use coaxial cables 
of time delay T coax in the control loop to approximately 
match the propagation time of the control signal with 
Tf. We denote this net control-loop delay time by r c = 
CLT + r coax (see Fig. [J) and adjust its value to be r c ~ 
40 ns (See Appendix for a detailed description of the total 
control loop CLT). 

In the following sections, we analyze the dynamics in 
v(t) from the full system (shown in Fig. [T]) and its rela- 
tion to the digital output of the control loop s(t). 



IV. FULL SYSTEM RESULTS 

We initialize the full system with v c t\ = wl and g(t) = 
9h to ensure oscillation growth. We monitor the tem- 
poral evolution of v(t) and s(t) using an 8-GHz-analog- 
bandwidth 40-GS/s oscilloscope (DSO80804A) and an- 
alyze separately the time series and spectral content of 
both signals. Along with our analyses, we explain the 
observed dynamics based on the system design. 



A. Experimental time series 

After intial transients, the system dynamics show oscil- 
lations in v(t) with a time- varying amplitude and digital- 
like switchings in s(t). As shown in Fig. H^,, the evolution 
of v(t) is aperiodic and its amplitude remains approxi- 
mately between ±50 mV, within the VGA's linear mode- 
of-opcration. The amplitude modulation in v(t) leads to 
spectral broadening of its power spectral density (PSD) 
with a peak at ~/ c and multi-mode power symmetrically 
distributed, as shown in Fig. 0t>. The small peaks in the 
PSD are spaced by ~l/rf. From the peak at f c , the PSD 
shows a 10-dB-bandwidth of ^100 MHz (consistent with 
the bandwidth of the LGA). Hence, the dynamics of v(t) 
have a UHF broadband spectrum. 

The complexity of the dynamics in v(t) has two differ- 
ent origins: (i) multi-mode feedback and (ii) long mem- 
ory. First, noise seeds the transient oscillations of the 
system, which are amplified and filtered based on the 
system's multi-mode characteristics (see the Appendix 



for the initial, transient dynamics). Then, these tran- 
sient amplitude modulations are stored in the system's 
memory as they propagate around the feedback loop and 
continuously the future dynamics at time t + Tf , causing 
the dynamics to reside in an infinite-dimensional phase 
space and allowing for the existence of chaos <21 

Simultaneously, the control loop generates s(t), a 
digital-like, asynchronous (unclocked) voltage. It 
switches apcriodically between its low and high states 
with minimum pulse widths of approximately 10 ns (see 
Fig. Ht). The finite bandwidth of s(t) is visible on the 
experimental PSD spanning frequencies from dc to ^100 
MHz, the cuttoff frequency at -10 dB (see Fig. Hi). The 
signal s(t) also shows short-pulse rejection^ a mecha- 
nism by which the TTL gate generates intermediate volt- 
ages in response to fast fluctuations at its input. As a re- 
sult, our system shows slight deviations from ideal PWL 
behavior, discussed in the following section. 

We finish this section by comparing the time series of 
v(t) and s(t). Unlike the audio-frequency PWL system 
presented in Ref . [l2| , there is no obvious correlation be- 
tween the analog and digital dynamics. This feature is 
unique to our system and is caused by the multi-mode 
feedback loop. In the next subsection, we construct a 
time-delayed return map of the analog dynamics to bet- 
ter understand the connection between v(t) and s(t) and 
reveal the symbolic nature of s(t). The return map also 
provides a method for calculating a Lyapunov exponent 
to confirm the presence of chaos. 



B. Experimental Return map 



In our system, shown in Fig. [TJ we note that 



>(*) = g(t)v(t - T f ). 



(2) 



We construct a return map by discretizing Eq. ^ . Sim- 
ilar to Ref. 12(, we use the local maxima in v(t), labeled 



as v n , where n is an integer incremented over each max- 
ima, that are spaced by ~l// c , as shown in Fig. [5^. Also 
shown is the local maxima i> n +T, where T is the integer 
number of local maxima that are contained in the feed- 
back loop. Using these points, Eq. becomes a delayed 
return map, defined as 



v n+T = M(v n ), 



(3) 



where M(x) = g n x and g n is the net gain of the system 
at time t n . 

The map is resprcscntcd graphically by a density plot 
in the plane of [v n ,v n +T)- Figure [SJd shows an attrac- 
tor of two different clusters of points with finite widths 
and slopes that coincide approximately with the gains gn 
and 5l (labeled accordingly in the figure). The linearity 
of these clusters is expected because the derivative of the 
map M'{x) = dM(x)/dx = g n . The finite widths relate 
to small over and under shoots in the values of the gain 
because of imperfections in the VGA gain and intermedi- 
ate s(t) values applied to v ct i- In addition, depending on 
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FIG. 4. Experimental time series and power spectral densities (PSD's). (a) Time series of v(t) showing UHF oscillations with 
complex amplitude modulations. A zoom of v(t) is shown for t = 0.5 — 0.525 /Lis. (b) PSD of v(t). (c) Time series and (d) 
PSD of the switching state s(t). The PSD's in (b) and (d) have a resolution bandwidth of 80 kHz and are smoothed over a 
window of 10 MHz. We note that the dynamics for this system are stable for extended periods of time, even when no thermal 
stabilization is used to prevent parameter drifts. 



the state of s(t), the map's dynamics reside on one of its 
two clusters, respectively. Everytime s(t) switches, the 
map's dynamical state (v n , iVi+t) transitions from one 
cluster to the other, demonstrating that s(t) partitions 
the map's attractor (See the Appendix for a preliminary 
analysis of the alphabet in the symbolic dynamics and 
its underlying grammar). 

We also analyze the map by looking at the statistics 
of its derivative. In Fig. [SJ;, the probability density 
of M'(x), labeled by Pm'(x), is given as a function of 
M'(x) and shows two primary peaks that occur at ap- 
proximately <7l,h, the most probable values of g n . There 
is also a higher probability to have M'(x) > 1 because 
of the larger area under the <?h peak, which indicates on 
average more growth than decay in the dynamics. 

Based on the finite switching times in s(t), we expect 
finite, but small, probabilities Pm'(x) ^ <?l,h, which is 
verified in the plot. However, the slow drop-off in Pm'(x) 
for M'(x) > <7h shows that there is a significant probabil- 
ity to have a gain g n > gn . This shows that the VGA has 
a slightly asymmetric response for switching v c t\ between 
Vh and i>h- Together, these features reveal the effects of 
the slight deviations from an exactly PWL system. 

To conlude our analysis, we use the probability density 
of M'(x) to calculate the Lyapunov exponent of the map 
(similar to Ref. (2(|) given by 

h = J\n\M'(x)\P M , [x) &x. (4) 

We note that, for ln|M'(a;)| > 0, M'(x) > 1 corresponds 
to exponential growth in the map's dynamics. Using 
Eq. @, we calculate h = (7.2 ± 0.2) x 10~ 3 , where the 
uncertainties are due to experimental noise and h > 
is an indication of chaos (using surrogate data of equal 



length, we have verified that h is significantly positive). 
In this case, h describes the average rate of divergence 
between any given pair (v n , v n+ x)- For the divergence 
rate between v n and a neighboring point u n +i, we scale 
h -> h/T with T = 87 such that h/T = (8.3±0.2) x 1(T 5 . 
With both Lyapunov exponents positive, this confirms 
the presence of chaos in our experimental system. 



V. PIECEWISE-LINEAR PHYSICAL MODEL 

In this section, we model both the feedback and control 
loops using using delay differential equations (DDE's). 
All of the model parameters are obtained from regres- 
sion analysis using experimental data as detailed in the 
Appendix. 



A. Multi-mode Feedback Loop Model 

We model the feedback loop as a band-pass filter with 
time-delayed feedback and gair>2£ 

v(t)+A(t)v(t)+^°\t) fv(t') dt' = g{t)v(t - r f ), (5) 

where rf is the feedback-loop time delay and A(t), uj(°'(t), 
and g(t) are approximately piecewise- constant (APC) pa- 
rameters that continuously switch. We note that Eq. ([5]) 
is completely linear for fixed A(t), o/°)(i), and g(t), as we 
have defined the feedback loop without a nonlinear sat- 
uration term. Therefore, for g(t) > 1, Eq. (|S|) becomes 
unstable and v(t) diverges. As we will discuss in sub- 
section B, the control loop, which continuously switches 
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FIG. 5. Construction of the experimental return map. (a) 
A small segment indicating the local maxima v„ (circles) to 
Vn+T, where T = 87. (b) Return map of v„. The dashed 
line indicates a reference line with slope equal to 1. (c) The 
probability density Pm'(x) as a function of AI'(x). The two 
primary peaks in Pm'(x) occur at M'(x) ~ <?l,h, the gain- 
states of the system. 



the values of the APC parameters, is the only mechanism 
that keeps v(t) from diverging (experimental saturation). 

The first two APC parameters are modeled as uj^ (i) = 
2n^/f(+)(t)f(-)(t) and A(f) = 27r(/<+>(t) - /<")(<)), 
where and /( _ ) are the upper (+) and lower (— ) 
cutoff frequencies of the multi-mode bandpass filter (-3 
dBm drop off), respectively. In the model, these cuttoff 
frequencies shift depending on the switching state s(t) as 



= ^ (ft''" 



- ft-') + fk + '~\ (6) 



where 



2.41 ± 0.01 GHz, ) 



1.85 ±0.01 GHz, 

1.97 ±0.04 GHz and 
s(t)/A is goverened by the model for the control loop. 
The last APC parameter is the gain g(t), which is also 



f[ +) = 3.29 ± 0.1 GHz, and f[ ) 



dependent on s(t)/A through the relation 
s(t) 



(7) 



where gn = 1.5 (<?l = 0.2) is the high (low) gain state. In 
our model, all of the APC constants have the same rise 
and fall times (instead of infinitely fast switches) that fol- 
low the physical switching state s(t). This approximation 
captures some of the non-ideal behaviors of the experi- 
ment while maintaining approximate PWL dynamics. 



B. Control Loop Model 

The control loop is modeled using two separate non- 
linear operators to approximate the LGA and TTL gate. 
The LGA first rectifies a delayed version of v(t)/4, pro- 
vided by the feedback loop, with a logarithmic response 
in amplitude and saturation for high inputs using the 
nonlinear function 



F(v(t-T c )) = Cx-C-i tanh 



log 



10 



-j-v{t - T c 



1 



(8) 

with parameters d = 1.68 ± 0.01 V, C 2 = 3.24 ± 0.09 V, 
C 3 = 24.1 ± 5.6 V -1 , and t c = 39.9 ± 1 ns. To finish, the 
output of this function drives a first order low-pass filter 



i(t) = -2nte(t) + 2nF{v{t - t-c)), 



(9) 



where /l = 0.02 ± 0.01 GHz is the low-pass cut-off fre- 
quency The output is smoothed to recover the en- 
velope of v(t). This completes the model of the LGA. 

The LGA output then drives the TTL gate, which 
we approximate as a continuous nonlinear switching func- 
tion 

s(*) = |(l+tanh[m(7 -£(*))]), (10) 

where A = 2.00 ± 0.04 V, I = 0.96 ± 0.01 V, and 
to = 51.00 ± 0.01 V -1 . Equation (HQ) switches with 
rise and fall times proportional to m and asymptotically 
approaches s(t) = and s(t) = A for inputs £(t) above 
and below the threshold /, respectively (Recall that the 
value of s(t)/A is also used in Eqs. (O - (0). 



C. Simulation Results 
1. Numerical time series 

We integrate Eqs. <[5|)-([7|) using a third-order Adams- 
Bashforth algorithm. The simulated time series and fre- 
quency spectra for v(t) and s(t) are plotted in Fig. [HI 
which mirrors the format of Fig. [4] for easy comparison 
with the experimental results. 

As shown in Figs. [5Ji-b, the central frequency of v(t) is 
~/ c , and the amplitude of v(t) is bounded between ±80 
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FIG. 6. Dynamics from the model, (a) Time series of v(t) with a zoom shown for t — 0.5 — 0.525 fis and (b) PSD of v(t). (c) 
Time series and (d) PSD of s(t). The integration was performed using initial conditions for t < : s(t) = 0, £(t) — 1.5 and 
v(t) — Gaussian white noise with a variance a = 10 -9 V 2 and sampled at a fixed time-step St — 5 ps. No additional noise is 
introduced at later times into the model. The PSD's are plotted assuming 50 Q impedance with a frequency resolution of 80 
kHz and averaged over a window of 10 MHz. 



mV. Similar to the experiment, the amplitude modula- 
tions of v(t) are aperiodic and the PSD of v(t) shows a 
multi-mode power spectrum. The dynamics of the model 
are thus a good approximation to the experiment, pro- 
ducing a qualitatively similar, broadband, UHF, multi- 
mode v(t). 

Recall that the model does not include any saturation, 
and thus, in order for such dynamics in v(t) to remain 
bounded, the control loop must produce s(t) to swtich 
the gain. In Fig. [l)p, s(t) is an asynchronous, digital- 
like waveform with some instances of short-pulse rejec- 
tion due to the continuous nature of Eq. (fTUj) . Different 
from the experiment, the PSD of s(t) shows small os- 
cillations in its amplitude. We hypothesize that these 
oscillations are due to high-frequencies in v(t) that bleed 
through the model's filter for the LGA. A higher order 
low-pass filter in the model can suppress these oscillations 
but requires an increase in the model's complexity. 

Next, we examine the PWL nature of the physical 
model by constructing a time-delayed return map of the 
integrated dynamics in v(t). 



2. Simulated return map 

Similar to the analysis of the experiment, we construct 
a time-delayed return map for v(t) using Eq. ([3]). Fig- 
ure [7^ shows that the dynamics of this return map form 
an attractor with two linear clusters that show qualita- 
tive agreement with the experiment, where the clusters 
produced by the model are less scattered. 

These quantitative differences occur because our model 
does not include noise and assumes the gain to be always 
ideally and linearly related to g(t). Nevertheless, in both 




40 80 

v n (mV) 

(b) xio- 3 




1 2 



FIG. 7. (a) Return map from simulation constructed from v n 
and v n +T, where T = 87. A dashed line indicates a reference 
line with slope equal to 1. (b) The probability density Pm'(x) 
for M'{x). 



the experiment and model, the return map transitions 
between the <?l and gn clusters as s(t) switches states. 

We also examine the simulated probability density 
-Pw'(x)) given in Fig. [7p. It shows two primary peaks 
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located at M'{x) ~ #l.h- Thus, the simulated dynam- 
ics also have high probabilities for the gain to satisfy 
g n = 9l.h and little probability for other values. Differ- 
ent from the experiment, these two peaks have smaller 
widths due to the assumptions of the model. 

Lastly, we use the simulated Pm'(x) and Eq. (j4j to 
calculate the map's Lyapunov exponent. We calculate 
h = 5.6 x 10" 3 , and the scaled h -> h/T = 6.4 x lO" 5 . 
We note that these values are close to the Lyapunov ex- 
ponents determined from the experimental data and con- 
firm the presence of chaos in the model. 

Taken together, the time series analyses and return 
map demonstrate that our simple model captures the 
fundamental aspects of the experimental system. 

VI. CONCLUSIONS AND FUTURE DIRECTIONS 

In conclusion, we demonstrate a UHF electronic de- 
vice displaying PWL chaos using feedback and control 
loops that overcomes many of the hindrances of slower 
LRC oscillators. We exploit a simple VGA with time de- 
layed feedback to generate growing and decaying multi- 
mode oscillations centered at frequency —2.1 GHz. Us- 
ing matched time delays for the feedback and control 
loops, our circuit performs fast control and avoids non- 
linear saturation. As for the dynamics, this novel cir- 
cuit generates simultaneously a chaotic signal with a 
corresponding switching state that partitions the return 
map's attractor, thus demonstrating a UHF foundation 
for the concepts of PWL chaos with a real-time symbolic 
dynamics . 12 i 19 

We presented our experimental findings alongside a 
numerical model. To capture the key aspects from the 
experimental dynamics, we measure the circuit's time- 
delays and frequency dependent band-pass characteris- 
tics. We simulate the dynamics using delay differential 
equations involving a logarithmic function, a low-pass fil- 
ter, and continuous switching functions to approximate 
the piecewise-constant parameters. Our model demon- 
strates a fundamental understanding of the system dy- 
namics and gives insight into the non-ideal effects arising 
in the experiment. 

In the future, we are interested in applying this fully- 
electronic design in a radar system. Our chaotic sys- 
tem automatically produces an upconvertcd, broadband 
analog waveform, and information about this signal can 
be stored with one-bit digital sampling of s(t). Prelim- 
inary analysis shows that the waveforms v(t) and s(t) 
show potential as radar signal sources when compared to 
standard techniques^ Futhermore, based on the PWL 
nature of this system, we are currently investigating a 
matchco™ (or pseudo-matched^) filter for chaos that 
can recover s(t) from v(t) in the presence of large noise. 
The simplicity of our design combined with these bene- 
fits represent a significant step towards inexpensive and 
robust chaos-based radars. 

Finally, we present only the simples case for our sys- 



tem where r/ and r c are approximately matched. Pre- 
liminary experimental and numerical work shows that 
additional mismatch between these two delays can give 
rise to new chaotic dynamics with differently structured 
PSD's. Thus, the delay mismatch represents an un- 
explored degrce-of-freedom for controlling the dynamics 
and exploring the stability of this system i 22 ' 30 Potential 
applications for this additional control include realizing a 
chaos communications syste m 31 ' 32 and creating orthogo- 
nal communications channels^ 3 - 
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Appendix A: Control Loop Characterization and Latency 
Time 

In this section, we characterize the three main com- 
ponents of the control loop: the LGA, the TTL gate, 
and the voltage control port v c t\ of the VGA, in order to 
estimate model parameters and approximate each com- 
ponent's control latency time (CLT). 

In Fig. [S^i, we plot a typical output waveform from the 
LGA. Using a signal generator, we drive the LGA with 
a 2 GHz sinusoidal signal x(t) that is pulsed for 10 ns. 
When zero power is present in x(t), the LGA outputs 
a positive voltage ~ 1.5 V. For oscillations in x(t), 
decreases to a value that is proportional to the am- 
plitude of x(t). Hence, the output of the LGA inversely 
follows the envelope of input oscillations. These time se- 
ries are measured using two high-impedance probes with 
identical propagation times after calibration. Thus, the 
time skew between x[t) and accurately represents the 
causal behavior of the LGA and shows that the CLTlga 
is approximately 8 ns. 

Using Eqs. ©, we drive our model for the LGA 
using x(t) and fit the output with the model pa- 
rameters Ci (for i — 1,2,3) and /l using a reduced x 2 
algorithm. We note that during the fitting, the CLT of 
our LGA model is not 8 ns due to the choice of a low- 
order filter. Thus, extra delay was added into the model 
LGA's CLT. This extra delay will be discussed at the end 
of this section. The resulting fit is shown in Fig. [SJd. 

We next characterize the TTL gate using the output 
from the LGA when it is driven by x(t). Using the same 
high-impedance probes, we measure and plot the TTL 
output s(t) alongside its input £(i) in Fig. \8t- In the 
figure, s(t) switches high (low) for £(t) below (above) the 
TTL threshold of -IV. We mark CLT T tl = 1-1 ns, the 
average delayed response of the gate switchings (rise and 
fall) in response to a threshold crossing. 

We use Eq. (fT0|) to fit the TTL gate parameters A, m, 
and /. The simple model does not produce the correct 
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FIG. 8. Control loop fits and latency times, (a) Output 
(top) of the LGA due to x(t) (bottom). The CLTlga 
is shown as a delay time between the two waveforms, (b) 
Model fit £(i) (top) due to input experimental x(t) (bottom), 
(c) Example output s(t) (green) of the TTL gate due to an 
input £(t) from the LGA. The TTL gate threshold is given 
by the horizontal bar. A typical CLTttl is shown for the 
switching in s(t) from high to low. (d) Model fit s(t) due to 
input experimental waveform £(t). 



CLT for the TTL gate, as an additional 1.1 ns time delay 
is necessary to achieve a good fit (see discussion at the 
end of this section). The fit is shown in Fig. [SJi. 

Finally, we analyze the processing time of the voltage 
control port of the VGA in reponse to a signal from the 
TTL gate. To do so, we drive the VGA with a 2 GHz 
low-amplitude sinusoidal signal and observe the output as 
s(t) switches to measure CLT 1 , ctl ~ 4 ns. Since the model 
for the gain of the VGA does not include a processing 
time, additional delay is necessary in the model. 

Therefore, in order to match and model the full control 
loop delay time r c , we use the relation t c = CLT + T coax , 
where CLT = CLTlga + CLT TT l + CLT„ ctl . We model 
the coaxial cable delay as a time shift that compensates 
for any CLT that is not present in our models for the 
LGA, TTL gate, or voltage control port. As a result, we 
tune the value of r c in the model to match that of the 
experiment, which is approximately 40 ns. 




f (GHz 



FIG. 9. Using Eq. (|B2[) . we fit the experimental spectra 
measured in Fig. [2] for (a) v c ti = and (b) « c ti = Wri- 



the setup shown in Fig. ®l using the differential equation 



m 

A l ,h 



-"(*)- 



A l .h 



= ah,Jl(v(t—T[)+Vin.), (Bl) 



where <zl,h = 9o,9l,yi is the attenuated gain (to keep the 
system from saturating) and Al,h, ^lh are the bandpass 
filter parameters for low (L) and high (H) gain states. We 
note that this equation is similar to that Eq. (J5J) with an 
additional input driving term Wj n and an attenuation g a 
in the feedback. We then Fourier transform Eq. (|BL[) and 
solve for £Tl,h(/) = ^/^im the system's transfer function, 
where v and Vi n are the Fourier transforms of v(t) and 
Vi n (t), respectively. The functional form of Hl,h(/) is 



#l,h(/) 



a], hAl.h 



2m f + A L , H 



aL.eA L . H e 211 /^ 
(B2) 



where Al.h = 2tt(/$ - /^), and = yj f^f££. 

Using the magnitude of i?L,H(/)j we fit the exper- 
imental data and obtain the fits shown in Figs. [9p,- 
b. The fitting parameters are /^ +) = 2.41 ± 0.01 GHz, 
= 1.85 ± 0.01 GHz, f[ +) = 3.29 ±0.1 GHz, and 
f[~^ = 1.97 ± 0.04 GHz. We note that the uncertainties 

are higher for the frequency parameters f^~' due to ap- 
proximation of the resonances in Fig. using Eq. (|B2[) . 
The experimental system is not a true band-pass filter for 
v c t\ = uh, and our approximation leads to discrepancies 
between Fig. f2t and Fig. [5p, but this approximation is 
necessary for our simple model. The fitted gain param- 
eters are an = 0.6 ± 0.01 and ol = 0.2 ± 0.01. Also, in 
order to facilitate fitting, tj is a free fitting parameter, 
where its fitted values are r/ = 41 ± 1 ns. 



Appendix B: Multi-mode feedback transfer function and 
fitting 

We determine the value of /l^ - " 1 (used in Eq. ©) 
using a fit to the experimental spectral responses shown 
in Figs. [5}>c. To do so, we derive the transfer function for 



Appendix C: Transient dynamics 

We plot typical transient dynamics for v(t) and s(t) in 
Fig. [TDJ Electrical and additive noise seed the experi- 
ment (Fig. [TUb ) and simulation (Fig. HUb). respectively, 
which show multi-mode growth in v(t) and that s(t) be- 
gins switching at 0.5 /is < t < 0.75 /is. 
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FIG. 10. Transients in v(t) and s(t) for the (a) experiment 
and (b) simulation. 



Appendix D: Underlying Grammar 

To demonstrate the existence of a potential grammar 
in the experimental and simulated switching states, we 
examine four-bit digital words that appear sequentially 
in the time series for s(t). To do so, we sample s(t) at a 
fixed clock frequency of 1 GHz and convert the digital- 
like signal into a list of digital ones and zeros, labaled s n : 
for s(t n ) > l,s n = 1 and for s(t n ) < 1, s n = 0, where 
t n are the clock sampling times. Next, we construct the 
pairings [0,0], [0,1], [1,0], and [1,1] by pairing [s n ,s n+1 ], 
[s„+2, s„+3], ... and then follow the evolution of the pairs. 
We depict the sequences of four-bit words in Fig. [TT1 
where a point is drawn between two pairings on the (x, y) 
axes if the four-bit word ([s„, s„+i], [s„ + 2, s n+ 3]) exists 
in the alphabet. Lines connect sequential four-bit words. 
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[0 1] 
[0 0] 
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[1 0] 
[0 1] 
[0 0] 



[0 0][0 1][1 0][1 1] 

X 



[0 0][0 1][1 0][1 1] 

X 



FIG. 11. Four-bit digital word maps for 10 4 words in the (a) 
experimental and (b) simulated s n . 

In the experimental analysis shown in Fig. II lb ,, there 
are both missing words and missing links, showing a re- 



stricted alphabet and grammar. In the analysis of the 
simulation (Fig. II lb ) . all four-bit words are present in 
the alphabet but with a restricted grammar. The gram- 
mar of the model is sensitive to parameter changes, and 
we conjecture that mismatches with the experimental pa- 
rameters, such as the order of the low-pass filter in the 
LGA, lead to the differences between Fig. [TTk and Fig. 
[TTb . We also note that overlapping lines are neglected 
and would show even more missing links in a higher- 
dimensional pictoral representation. 
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